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ABSTRACT 

Position and frequency switching techniques used for the removal of the bandpass dependence of radio astronomical spectra are 
presented and discussed in detail. Both methods are widely used, although the frequency dependence of the system temperature 
and/or noise diode is often neglected. This leads to systematic errors in the calibration that potentially have a significant impact on 
scientific results, especially when using large-bandwidth receivers or performing statistical analyses. We present methods to derive 
an unbiased calibration using a noise diode, which is part of many heterodyne receivers. We compare the proposed methods and 
describe the advantages and bottlenecks of the various approaches. Monte Carlo simulations are used to qualitatively investigate both 
systematics and the error distribution of the reconstructed flux estimates about the correct flux values for the new methods but also 
the 'classical' case. Finally, the determination of the frequency-dependent noise temperature of the calibration diode using hot-cold 
measurements or observations of well-known continuum sources is also briefly discussed. 
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1. Introduction 
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Spectroscopic data obtained from radio astronomical observa- 
tions typically need to be post-processed in several steps be- 
fore one can proceed with the scientific analysis. The uncali- 
brated flux values not only have to be converted from device 
units (counts) to physical quantities (e.g. Kelvin or Jy/beam), 
but one must also deal with the frequency-dependent gain (band- 
pass curve), radio frequency interference, removal of any resid- 
ual baselines, and so on. While several me thods are disc ussed in 
the literature to solve the bandpass issue (lHeilesll2Q07l and ref- 
erences therein), the flux calibration is often handled using only 
a very simplified approach. In particular, the widely used posi- 
tion and frequency switching techniques (see Heiles 2007, for a 
review) suflTer from an improper treatment of calibration. In this 
paper we describe the underlying problem, present possible so- 
lutions, and assess the quality of the various approaches in terms 
of flux calibration accuracy. We compare the proposed methods 
with the 'standard' approach, and show that the latter is subject 
to bias eflTects that can have a large impact on scientific results. 

The paper is organised as follows. Section|2] introduces the 
basic equations and definitions necessary to describe the prob- 
lem. Sections[3]and|4lpresent the position and frequency switch- 
ing procedures, along with several methods to correctly cali- 
brate the measured spectra. An alternative approach is shown 
in Section[5] which may be of interest in cases where neither of 
the switching schemes are applicable. To illustrate the calibra- 
tion methods, artificial spectra were generated with well-defined 
input quantities. Section[6] discusses the various approaches and 
their ability to reconstruct the simulated flux values. The flux cal- 
ibration presented here relies on the signal of a noise diode fed 
into the receiver. Two possibilities for measuring the spectrum 
of this calibration signal are discussed in SectionQ We conclude 
with a summary in Section[8l 



2. Basic equations and definitions 

The output signal, P, of a spectroscopic (heterodyne) receiving 
system can be described with the simple formula 



P^'^Hv) = GRF(y)GiF(y) (rsou(v) + ^^^^"(v)) , 



(1) 



where the telescope itself and the receiving system exhibit a 
frequency-dependent gain, G = GrfGif, on the incident radi- 
ation. Often referred to as Ta, T^qu, is the antenna temperature 
caused by the flux density of the astronomical source of inter- 
est, which potentially incorporates line emission and a contin- 
uum contribution, i.e. T^ou = T^^Jfu + T^ou- Apart from the flux of 
the astronomical source, Tsou, that we are interested in, there are 
several noise contributions, which we subsume in the so-called 
system temperature, Tlyf^ Q This consists of 

^sys - ^bg + T'atm + T'spiH + T^sw + T'loss + T'rxL+T'cal] (2) 

Tbg . . . microwave and galactic backgrounds 
Tatm . . . atmospheric emission 
Tspiw . . . ground radiation (spillover and scattering) 

Tsw . . . standing wave pattern (often in secondary focus) 
Tioss . . . losses in feed, ohmic losses 

Trx . . . receiver noise temperature 

Teal . . . injected noise using a noise tube/diode. 

Standing waves (SW) are produced when the incident radio sig- 
nal is reflected at the aperture of the antenna leading to in- 
terference patterns that often exhibit a sinusoidal sha pe. For a 
more detailed discussion, we refer the reader to iRohlfs & Wilsonl 
(120041) . 



^ Note that sometimes the system temperature is defined diff'erently 
to include Tsou- 
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In this paper we make extensive use of the calibration sig- 
nal Teal, which is a part of many heterodyne radio receivers. Its 
purpose is to provide a well-defined intensity standard that can 
be used to calibrate the measured spectra in terms of the antenna 
temperature. In practice, there are two common techniques to 
make use of Tcai- One is to rapidly switch the Tcai signal on 
and off' on time- scales of the order of tens of milliseconds up 
to a few seconds throughout the observation, leading to two so- 
called switching phases cat and non-cal. This method is used at 
the 100-m Eff'elsberg telescope, for example. In the following, 
we adopt the notation P and P^^^ when referring to non-cal and 
cal phases respectively, whilst P^^^^^ indicates either of the two. 
The alternative method in use slowly cycles between the two 
phases during a dedicated measurement of a calibration source 
(or blank sky) every few hundred seconds, e.g. at the Arecibo 
telescope. In this paper, we assume a rapid Tcai switching for 
the sake of a simpler presentation of the calibration methods. 
However, our work may easily be extended to the slow cycling 
scheme. We would also like to point out that many (sub-)mm- 
wave telescopes make use of chopper wheels or hot-cold loads 
to improve the calibration. It should again be straightforward to 
adapt our methods to these techniques, though this is beyond the 
scope of this paper. 

A spectroscopic back-end usually records the spectral den- 
sity function, P, in arbitrary units (e.g. counts) such that the raw 
spectra need to be calibrated in terms of both ffux density val- 
ues and also with respect to the frequency-dependent bandpass 
shape G. We note that for simplification we incorporate the co- 
efficient that converts the antenna temperature to counts into the 
gain factor G. Using a heterodyne receiver, one has to distinguish 
between the gain in the receiver/radio frequency (RF) and in the 
intermediate frequency (IF) part of the receiving system. 

One major problem in the reduction of spectroscopic data 
is that the bandpass curve(s), G, must be disentangled from the 
input signal, T. From Eq. ([T]), it is clear that one has to divide 
the measured signal, P, by the bandpass, G, to obtain the in- 
put temperature spectra. Unfortunately, the bandpass curve is 
unknowr0. To cope with this issue, one often applies a switch- 
ing technique to obtain a reference spectrum, e.g. position or 
frequency switching as explained below. Dividing by this refer- 
ence removes the frequency-dependent gain. Residual baselines 
caused by SW, for example, may then be subtracted. We note 
that one can neither bandpass-calibrate a measurement by sim- 
ply subtracting a reference spectrum nor remove residual base- 
lines by division of a baseline model. 

Even after the (successful) removal of the gain curve, the 
extracted Tsou spectrum is not equal to the true brightness tem- 
perature, r*Q^, of the astronomical object. Both the atmospheric 
dampening and aperture efficiency, //ap, which act on r*^^ prior 
to amplification, must also be taken into account 



T/apr^ou^xpC-'^oAM), 



(3) 



where tq is the zenith opacity and AM = sin" ^ (El) the elevation- 
dependent airmass.0The aperture efficiency, T/ap, is the product 
of several other efficiencies. The antenna pattern and blockage 
of the aperture play a role, as well as the surface accuracy and 



^ Even if it had been measured, it is still subject to instabilities, e.g. 
caused by temperature drifts, such that frequent re-evaluation would be 
needed. 

^ Note that sin" ^ (El) does not accurately describe the airmass 
at low elevation angles and must be replaced by a site-dependent 
model, including vertical va lues of density and refraction index (e.g., 
iMaddalena & Johnson|[2005h . 



to some extent ohmic losses. This cannot be calculated theoret- 
ically, especially for the large single-dish telescopes. It is pos- 
sible to model opacity (see Kraus 2009). In the best case, an 
atmospheric model makes use of simultaneous measurements 
of water vapour to achieve greater accuracy (e.g. using a water 
vapour radiometer, iRottmann & Royll2007l) . For the remainder 
of this paper, we neglect these corrections for simplification as 
they are completely independent of the calibration methods pre- 
sented here. 



3. Position switching 

In position switching, one points the telescope first to a refer- 
ence position (Off) before measuring toward the true target (On). 
Applying Eq. ([T]), we obtain 

p[-" = GrfGif^-" = GrfGif (r,ou + r[-;]„) , (4) 



r)[cal] _ 
off 



n[cal] _ 



n[cal] 






(5) 



The aim of position switching is to remove the bandpass curve. 
In the ideal case, GrfGif is equal for both positions and dividing 
both spectra leads to 



p[cal] J . 7^ [cal] 

/^nn ^ SOU ^ ^ SVS.O 



p[cal] 

^Off 

This can also be written as 

p[cal] _ p[cal] 
^on ^off 



sou ^ ^ sys,on 
p[cal] 
" sys,off 



r.™, + Ar, 



sys 



p[cal] 
^off 



T-'[cal] 
sys, off 



(6) 



(7) 



where we introduce as 

Ar^ys = T^t, - T^^^^ ^ ATbg + Ar^ky + ATspii, + Ar^w, (8) 

the diff'erence between the system temperatures in the On and 
Off positions. The receiving system is usually stable for the 
duration of a single position switch and only a few contribu- 
tors must be taken into account. While Tsky and Tspm are ele- 
vation dependent, T^^ may be a function of the incident contin- 
uum flux. In some cases, SW would only be noticeable if one 
pointed toward strong continuum sources, e.g. calibrators, and 
in other cases during daytime because of the sun, regardless of 
the observed source. The Galactic background continuum radi- 
ation depends slightly on the pointing position of the telescope. 
Using an appropriate Off position, for instance at constant ele- 
vation angle, one should usually be able to minimise AJsys. The 
Off position should obviously not contain a continuum source it- 
self. The equations could otherwise still be interpreted in a way 
where Tl^^ is not the true continuum flux of the On position. 
Re-arranging Eq. (|7]), we obtain 



D[cal] 



-t sou ~r ^-t svs — -t , 



sys 



[cal] '_ 
sys, off 



-P' 



[cal] 
off 



D[cal] 
off 



or alternatively 



+ r 



[cal] 
sys, off 



+ AJsys = t 



_ ^[cal] 



3[cal] 



" sys,off p[cal] ' 
^off 



(9) 



(10) 



It is important to realise that T^^^^ ^^ is a frequency -dependent 
function. One can ensure a proper flux calibration over the com- 
plete spectrum only if this is determined correctly. In the re- 
maining part of this section, we develop two methods to de- 
termine ^f foff(v). The first is easier to apply and more robust. 
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The second may be of interest under particular circumstances 
when using frequency switching and is also discussed in the 
case of position switching for completeness. Both approaches 
make use of the signal of the calibration diode (or noise tube). 
Furthermore, we discuss the 'classical' method, which assumes 
that Tsys = constant and Tcai = constant, as a simplification of 
the first method. 



3.1. Setting up simulations 

In the following subsections, we present the three alternative cal- 
ibration methods. Artificial spectra were generated to enable a 
direct comparison of input and output spectral line intensities. In 
Section[6l we statistically assess the quality of these methods to 
evaluate their general robustness and error distribution. 

The test spectrum uses a relatively simple set-up. The con- 
tinuum flux of the Off position is described as a single-slope 
power law rsys,off = 400(y/yo)"^-^ [K], where vq = 300 MHz. 
This is a strong simplification of reality where several diflTer- 
ent power laws (e.g. from the galactic background, sky, and 
ground radiation) add to a more or less 'arbitrary' tempera- 
ture contribution from the receiver. For the astronomical source, 
Tsou, a set of three Gaussians superposed on the power law 
r,^o7 = 200(y/yo)-^-^ [K] (where vq = 300 MHz) was used. 
Each Gaussian has equivalent line parameters (3 K amplitude, 
1.4 MHz line width, FWHM) but are added at diflTerent fre- 
quencies. These are used as 'probes' to measure the flux cali- 
bration quality at diflTerent frequencies after data reduction. We 
chose to model a possible frequency dependence of the noise 
diode using a shallow power law Tcai = 3(y/yo)"^'^ [K] (where 
vo = 1420 MHz). In reality T^ai is usually a much more com- 
plicated function especially at higher frequencies, but this setup 
is suitable for the main purpose of our simulations, which is to 
illustrate the various steps in the data processing as simply as 
possible. Section [T6l illustrates a more realistic test case with re- 
spect to the functional form of the system temperature and cal- 
ibration diode. The resulting temperature spectra are plotted in 
Fig. [T] (upper panel). 

Gaussian noise was added to each component according to 
the root mean sq uare (RMS) noise calcu lated using the radiome- 
ter equation fsee lRohlfs & Wilsonll2QQ4h 



Trr.. = AT 



kT 



input 



(11) 



where ^ is a constant factor depending on, e.g. the input quanti- 
sation of the signals. A/ is the bandwidth of a spectral channel, 
and T is the integration time. In the simulations, we assume the 
total number of spectral channels to be 16k over a bandwidth of 
300 MHz and an integration time (per spectral dump or measure- 
ment phase) of 5^. For example, a system temperature of 25 K 
would result in Tnns ~ 83 mK (per spectral channel). Eq. (fTTI) 
implies that noise is a function of frequency because the tem- 
perature values, Tinput, for each measurement phase vary over 
the observed band. Likewise, the additional temperature compo- 
nents of the continuum emission of the observed source and the 
calibration diode increase the noise. It is also possible to cal- 
culate the theoretical RMS level expected in the final (reduced) 
spectrum using the noise values of the individual input spectra 
(see AppendixlAb. 

Finally, for a more complicated test case, we simulated 
(mono-modal) SW by adding a sine wave component T^^ = 
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Fig. 1. To test the influence of the diflTerent calibration schemes 
on the reconstructed flux values, synthetic spectra were pro- 
duced. The upper panel shows the input power spectral densi- 
ties. The antenna temperature on-source (Tsou) is a superposition 
of three Gaussians with an amplitude of 3 K on top of a contin- 
uum spectrum with a spectral index of -0.7, i.e. T^ou ~ y~^'^. For 
the system temperature, rsys,off ~ v"^-^ is used for convenience. 
We also simulate a slight frequency-dependence of the calibra- 
tion diode using the power law Tcai ~ y"^'^. Gaussian noise 
was added to the four measurement phases according to their to- 
tal temperature values (see text). Note that each synthetic spec- 
trum contains a diflTerent noise realisation. In addition, a more 
complicated case is included with a standing wave contribution, 
which is parametrised as a monochromatic sine wave component 
(lower panel). 



2 sin(2;r(y/ysw)) [K] to rsys,off (see Fig.[Tl lower panel)o For the 
position switching test, ygw = 50 MHz. 

The input temperature spectra were multiplied with a sim- 
ple bandpass curve shown in Fig.O The resulting power spectral 
densities (in arbitrary units) are plotted in Fig. [3] The superpo- 
sition of the SW and the ripples in the bandpass shown in the 
lower panel of Fig. [3] creates a rather complex pattern. 



^ Note that this set-up does not account for all possible situations. 
For example, if strong continuum sources alone produce SW, then they 
contribute to 7^^°^^ only. However, this is only a minor problem as the 
computation of Tsys is not affected. Nevertheless, one has to account for 
such a case if one uses continuum sources for calibration. 



B. Winkel et al.: Unbiased flux calibration methods for spectral-line radio observations 




1300 1350 1400 1450 1500 1550 

V [MHz] 

Fig. 2. The bandpass curve, G = GrfGif, used for the simula- 
tions. 
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Fig. 3. Multiplying the input spectral densities with the bandpass 
shape shown in Fig. [T] returns the 'measured' spectra P. The up- 
per panel shows the results for the simpler case, while the lower 
panel includes the standing wave contribution. 



3.2. Using the Off position to obtain 


^ sys 


Computing 










peal 
off _ 


^ sys,off + ^ cal 


Teal 



/'off 



+ 1 



' sys, off 



- sys, off 



(12) 



leads to an equation to infer the system temperature of the Off 
position, which depends only on the temperature of the noise 
diode, Teal. If the latter is known, then 



- sys, off 



= T, 



cal 



peal 
£_off 

Poff 



-1 



(13) 



However, the noise diode temperature is a frequency-dependent 
quantity. It can be measured with good precision using the 
hot-cold method (compare also Section lTJl and Appendix IB.lK 
which, unfortunately, is a time-consuming procedure. One so- 
lution is to establish a catalogue of astronomical calibrators, 
i.e. bright continuum sources, which serve as reference to 
(re-)calibrate the Tcai spectrum on appropriate timescales (see 
Section lT!2l and Appendix lB.2b . This approach was also pro- 
posed by Maddalena & Johnson (2005). Since Tcai is a time- 
dependent quantity (e.g. owing to a change in the environmental 
conditions), it is always a potential source of (systematic) error. 
Nevertheless, a temporal stability on the order of 1% can typi- 
cally be expected on the timescale of one hour. 

Unfortunately, even if one has a good model of Tcai, the mea- 
sured quantity P^^/Poff is still subject to noise and will substan- 
tially increase the noise in the final reduced spectra. One solution 
is to suppress the noise in the obtained rsys,off spectrum before 
substituting it into Eq. Q. 

In simpler cases, namely in the absence of a standing wave 
contribution and a relatively flat rcai(v), the quantity 



peal 
f_off 



f sys, off 
^eal 



= Koff(y)- 



(14) 



is approximately proportional to rsys,off , which can be described 
by a power law. Hence, a low-order polynomial might already 
suffice to describe the quantity /Coff (v). 

Unfortunately, at higher frequencies in particular, it is almost 
impossible to engineer a sufficiently flat Tca^v). In such cases, 
a filtering approach or the use of high-order polynomials is re- 
quired to suppress noise (see Section [T6l for a more realistic ex- 
ample). In most cases, it is easier to model k~^ for numerical 
stability, as usually rsys,off » T^ai- The noise-free (or low-noise) 
model can then safely be inverted to obtain Kq^. 

Fig. 131 illustrates K~^(y) = 7'cai/7'sys,off for the two test cases 
with and without SW. In the former case (upper panel), a third - 
order polynomial was used to describe k~^. For the latter, a more 
complicated model is necessary 



^.^fF — 



Liinuv) 



'^off 



^2(^2, v) -h A sin(ay -h b) 



(15) 



where L/(^/, y) are polynomial functions of degree rit = 3. This 
result was obtained only after providing suitable initial fit pa- 
rameter values and is shown in the lower panel. An entirely auto- 
mated procedure would most likely have difficulties in handling 
SW — even for the simplest scenario of monochromatic SW. 
Substituting Eq. ([T4l) into Eq. ([9]) leads to 



^ sou"!"^^ sys ~ ^off ^ ( 



^on ^( 



off 



Deal 



f eal- 



PoS 



= (Koff+l)T, 



eal" 



peal 
^off 



peal 
^off 



(16) 



We note that the continuum flux of the source is also recon- 
structed. To obtain an optimal signal-to-noise ratio (S/N), the 
reduced cal and non-cal spectra should be averaged. 

After applying Eq. ([T6l) to the simulated data (using a poly- 
nomial model for /Coff; see Fig. (4]), the spectra shown in Fig. [5] 
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Fig. 4. Spectrum of k~^ for the simpler (upper panel) and stand- 
ing wave (lower panel) cases and the fitting models inferred. 



are obtained. To measure the spectral line intensities, we fit a 
Gaussian superposed on a third-order polynomial to each spec- 
tral line. We note that Fig. [5] shows the reduced spectra for the 
SW case. The non-SW scenario produces very similar results. 
The upper panel shows the direct result and the lower panel con- 
tains the spectrum after baseline removal. The noise values and 
the recovered line intensities match the expected results. The 
residual noise is higher than one might naively expect (the RMS 
in the final spectrum is around the same level as for each of the 
four contributing phases). This is due to the division by the ref- 
erence spectru m tha t does not contain signal but adds noise (see 
also Appendix lA.il) . 

It is possible to smooth the reference spectrum prior to divi- 
sion greatly reducing the RMS in the final spectrum by almost a 
factor of V2. However, as iBraatzl (l2009l) points out, this proce- 
dure can have unwanted side-eff'ects, for example, degradation 
of the baseline or emphasis of narrow features present in the 
reference, e.g. radio frequency interference (RFI), and results in 
(locally) correlated noise. 

One drawback of the proposed method is that the Off po- 
sition alone is utilised to infer T^y^. As a consequence, a factor 

of V2 in sensitivity is lost in calculating Koq. However, the final 
result is unaff'ected so long as Kqs is described by a (smooth) 
model. If sensitivity is not a major concern, even the origi- 
nal (noisy) Kqs spectrum, i.e., without modelling, may be used. 
This is advantageous since the solution is then unaff'ected by the 
choice of a given model. 
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Fig. 5. Upper panel: Resulting spectrum for the SW case after 
applying Eq. ([T6b , using the model fit shown in Fig.|4j The con- 
tinuum contribution of the source is recovered. A baseline is sub- 
tracted in the lower panel and spectral line intensities were fitted 
locally. The recovered intensities match the input very well. The 
noise values in diff'erent spectral regions are also consistent. 

3.3. Using On and Off positions to obtain T^ys 

The second method utilises Eq. d?]) and uses continuum contri- 
butions only, to define the quantities 



f(y)^ 



off 



r\y) = 



^ off 

peal _ peal 
^on ^off 



peal 
^off 



^ eont ^ ^^ sys 
^ sys, off 

^ eont ^ ^^ sys 



' sys, off 



+ 7-, 



(17) 
(18) 



eal 



Since the On spectra are involved, it is important that the fitting 
algorithm excludes spectral lines. In simple cases, the use of rel- 
atively low-order polynomial functions should suffice to provide 
a good model. If SW occur and only few modes (i.e. SW fre- 
quencies) are present, one should be able to model these with 
superposed sine-waves. In principle, one could even try to use a 
filtering approach to obtain a noise-free model /. 

The fits calculated for the example spectrum are shown in 
Fig.[6l For the non-SW case, a third-order polynomial model is 
used. The SW case uses 



f(y)^ 



Liinuv) 



^2(^2, v) + A sm(ay + b) 
where L/(^/, y) are polynomial functions of degree rii. 



(19) 
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Fig. 6. To correctly reconstruct the original fluxes, the sec- 
ond proposed method uses models of the intermediate spectra 
(Pof^ - ^lff")/^lff" (see Eq.dnl) and ^). Appropriate win- 
dows around the spectral lines should be set because one is 
only interested in the continuum contribution. The upper panel 
shows the results of the simpler case, utilising second-order 
polynomials sufl&cient to describe the baseline. For the stand- 
ing wave case (lower panel), a more complicated fitting model, 
Li(ni,y)/(L2(n2, v) + A sin(ay + /?)), was applied, where L/(f7/, y) 
are polynomial functions of degree rit. In this example, we use 
ni = n2 = 3. 



Substituting Eq. (fTTl) into Eq. ([TSl) gives 

/cal(y) 



^ sys,off ~ ^ ( 
^ sys,off "I" ^ cal — ^ cal 



"'/(v)-/^^'(v)' 

f(y) 
/(v)-/'=^'(v)' 



(20) 
(21) 



such that 7'sys_oif(v) is a function of Tcai. Eq. ^ now becomes 



yline , T^cont _i_ ax — T 
^ sou "•" '' sou "•" "^ sys — ■'car 



/cal(v) 



^on ^( 



off 



/(y)-/cal(y) p^^ 



Teal 



f(y) 



Dcal 



peal 
^off 



(22) 



/(y)-/cal(y) pc.1 



Using the calculated model fits and applying Eq. (l22l) , one 
obtains the reduced spectrum. Fig. [7] shows the standing wave 
example where the upper panel shows the reconstructed spec- 
trum, Tsou, and the lower panel contains the result after baseline 
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Fig. 7. Inserting the model /^^^^^ (v) into Eq. (1221) results in the 
correct flux calibration (upper panel). The lower panel shows 
the spectrum after baseline subtraction. Note that the correct 
continuum contribution of the source was also reconstructed (up- 
per panel). The example was computed for the SW case. 



subtraction. The continuum contribution, T^q^\ of the source is 
correctly reproduced. In some cases however, the baseline may 
contain a residual imprint caused by standing waves, which af- 
fects the quality of the model fitting to some extent, especially 
in the presence of noise. In the worst case, this might also af- 
fect the flux calibration. The more complex the required model, 
e.g. in the presence of SW, the larger the potential error in flux 
calibration. 

We note that this method may not be applicable in all cases. 
Inspecting Eq. ([22b in detail reveals that the denominator of the 
term 

/cal(y) 



/(y)-/cal(y) 



(23) 



can become zero. This occurs if T^^^^ -\- rsys,on 



^ sys, off? 1-e. 

for sources with low continuum flux. In these cases, a nearby 
calibrator may be used for an independent reference position in 
order to obtain a well-behaved correction term. 



3.4. Determination of Tsys negiecting any dependence on 
frequency — tlie 'classical' approach 

The 'classical' case where Tsys and Tcai are both treated as con- 
stants is clearly a simplified case of the first method (presented 
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in Section [T2l) and of course causes errors in the flux calibration. 
This simplification was justifiable to some extent until about two 
decades ago when broadband receiving systems started to be- 
come more commonly used. Today, this approach is still widely 
used, e.g. within the onl ine reduction pipeline at the 100-m tele- 
scope in EflTelsberg (see lKrausll2009l) . in order to provide a fast 
and robust online display of measured spectra. 
Eq. ([T2I) may be rewritten as 



- sys,off 
Teal 



^off 



(p--UP,s)-(P^-'-P,ff) 



off 



(^off-^off) 



2(/^off-^off) 



(24) 



In principle, this provides a direct means of calculating rsys,off 
from Teal. Unfortunately, a simple evaluation of the above equa- 
tion is numerically unstable because the denominator can have 
values close to zero or even become negative owing to noise. If 
we treat rsys,off and Tcai as a constant with respect to time and 
frequency, we find that 



' sys,off 
Teal 



(P-i + Poff) - (P'i - Po^}y 



2{Pli-Pos}r 



(25) 



where {P^^^ - P)y is obtained by calculating the average of the 
inner 50% of the diflTerence spectrum in the online reduction 
pipeline at the 100-m telescope at EflTel sberg. The c urrent data 
reduction pipeline at the GBT (GBTIDL; Braatz "200913 uses the 
equivalent method except using the inner 80% of the spectra. For 
future versions of GBTIDL, a vectorised approach is planned 
whereby averaging is performed on certain bins. 

The parameter Tsys is calculated by the same averaging pro- 
cedure, such that 



^ sys,off — V sys,off /y 



= r< 



call — 



2{Pl'^ - Poff)v 



(26) 



Using a scalar value for Tcai, roughly known for each receiver 
from Eq. (|25]) , we obtain the 'calibrated' spectrum 



■ ^ sys,off 



off /TF 



Dcal 



Poff 



{f 



sys,off 



+ rc 



cal 



peal 
■^off 



p^i 

off 



(27) 



Both, cal and non-cal spectra, should again be added to reduce 
noise in the final spectrum. 

Although this approach is much easier to implement in soft- 
ware, the reader should be warned that neglecting the frequency 
dependence of Tgys introduces systematic errors, which can be 
a serious problem. This is especially the case when computing 
line ratios or dealing with statistical analyses of large samples of 
sources at various radial velocities/frequencies (see Section[6]). 

Figure [8] (upper panels) shows the resulting spectra after ap- 
plying Eq. (1271) . We note that the baseline does not match the in- 
put baseline at all. For reference, we subtracted a baseline (which 
proved to be rather complicated in the SW case) from the result- 
ing spectra. The result clearly shows a systematic flux calibration 
error and is displayed in Fig. [8] (lower panels). 

3.5. Computing the bandpass curve G = GifGrf 
From Eq. (|5]), it directly follows that 

Poff 



G = GifGrf = 



peal 
^off 



- sys,off 



(^off + im 



(28) 



eal 
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where rsys,off can be calculated, for example, via Eq. ([T4l). 
7'sys,off = KoffTcah ov using Eq. (|2Q1) and (I2TI) . This can be of great 
practical use. While the inferred G from a single measurement 
(see Fig.O is rather noisy, one may use several observations, 
even of different sources, to average the gain curve and suppress 
noise. The calculated G could then in turn be applied to all in- 
volved data sets. In practice, the gain is usually imperfectly sta- 
ble in time. Eq. ([28l) enables us to monitor the time evolution of 
the receiving system. It is even likely that temporal drifts can 
be approximated to some degree with a low-order polynomial 
(spectral channel- wise) and can then be used to improve the ac- 
curacy of the data reduction. 

3.6. A realistic case 

We present the results of a more complex simulation that resem- 
bles realistic observations. In comparison to the previous exam- 
ple, both Tsys and Tcai now have more structure. Both quantities 
were modelled to imitate the outcome of a measurement with 
the 1.9-cm primary-focus receiver at the 100-m telescope. The 
resulting spectrum is illustrated in Fig.[TO] T^ys shows a steep in- 
crease towards higher frequency, attributed to the water vapour 
line at 22 GHz (during the observations, weather conditions were 
rather bad). 

In this instance, the first method as described in Section [T2l 
is applied using spectral filtering instead of polynomial models 
to describe a:^^. The result is shown in Fig.[TT] The second ap- 
proach requires spectral windows to be set around the (expected) 
emission lines such that a simple filtering approach would not 
work. It is therefore neglected here. 

Fig.[T2l displays the resulting spectrum. Despite the very 
complex set-up, the calibration worked well and the underlying 
continuum flux of the observed source could be reconstructed. 
Fig.[T3l shows the result of the classic approach. For purposes of 
comparison, a constant continuum contribution was chosen. The 
resulting baseline is very complex and is in fact the incorrectly 
calibrated continuum flux of the source. In a worst case scenario, 
this can lead to a misinterpretation of the results depending on 
the line widths of the expected astronomical features. For ex- 
ample. Fig. [13] suggests that the middle spectral line has a blue- 
shifted wing. 



4. Frequency switching 

In frequency switching, one uses two diff'erent local oscillator 
(LO) frequencies to provide a reference spectrum. The shifted 
spectrum serves as a reference to remove the bandpass depen- 
dence in the non- shifted spectrum. In practice, a symmetric LO 
shifting pattern is often used, such that Eq. ([T]) may be written as 

/^i-;"(v)=GiF(y)GRF(y-Ay)x 

[^sou(v-Ay) + r:-yy-Ay)], (29) 

/^5"(v)=GiF(y)GRF(y + Ay)x 



[rsou(y + Ay) + Tllfl^iv + Ay)] . (30) 



Together with the cal and non-cal phsiscs, the LO switching leads 
to a total of four diflferent phases of sig+cal, ref, sig, and refrcal, 
which are usually rapidly stepped through in order to avoid time- 
dependent instabilities. 

Using a shorter notation and assuming T^^^ \ (y) = T^^^ Jv), 

^ ^ sys,sig^ "^ sys,ref^ ^' 

which should be fulfilled since switching is performed on short 
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Fig. 8. Applying the 'classical' calibration scheme that is currently in use at the 100-m telescope (for online display purposes), using 
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We first assume that Grf,- = Grf,+, which is fulfilled only in the 
very rare cases that Grf(v) = constant. Then 



p[cal] _ p[cal] 
sig ref 

p[cal] 
^ref 



+ r 



[cal] 
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^ sys,+ 



+ r. 
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(33) 
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line,- 



I 'T-'SOU 
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line,- 
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where Ar. 



[cal] 

sys,i 



= r. 



[cal] 

sys,- 



- Tl^f\ = Arsys,± [+Arcai,±]. We note that 
Arsys,± is equal in the cal and non-cal phases, while Arcai,± is 
diflferent (Arcai,± = in the non-cal phase and Arcai,± 9^ in the 
cal phase). Both quantities depend only on the slope of Tgys and 
Teal because all contributors to Tsys remain approximately equal 
during two adjacent switching phases. Using Eq. (1331) , it follows 
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Fig. 10. Temperature inputs for a more realistic simulation. The Fig. 12. The resulting spectrum in comparison with the input 
three Gaussian emission lines have a 1 K amplitude and a line flux, 
width of 240 MHz (FWHM). The number of spectral channels is 
128k. 
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p[cal] _ p[cal] 
[cal] \ sig ref 



r,[cal] 
ref 



_ y^[cal] 
~ sig ' 



(34) 
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In Sections l43l and 1441 we discuss how to perform the calibra- 
tion for frequency switching, i.e. how to determine r^^ou,* + ^sys!±- 
To obtain an optimal S/N, one should not only average cal 
and non-cal phases but also shift (to match RF frequencies) and 



average f sig and f ^f 

lff5"(v + Av) + fr(v-Av)l 



T'souCv) - -TsouCv + 2Av) - -rsou(v - 2Av) (36) 



+ T^f\y) - \T'^f\y + 2Av) - l7-[-'I(v - 2Av). 



sys 



This is clearly a much more complicated procedure than for the 
position switching case. 

4.1. Bandpass ghosts 

In Eq. ([36b, apart from the desired signal, T^o^iv), so-called 
'spectral-line ghosts' (sometimes also referred to as bandpass 
ghosts), -jTsou(y ± 2 Ay), appear in the resulting equation. They 
have half of the amplitude of the original signal. Furthermore, 
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Fig. 14. Applying the frequency- switching base equation, 
Eq. (l36l) , one would expect the spectral-line ghosts to have half 
the amplitude of the original line. However, in practice one is 



.+r. 



[cal] 



rline 



n[cal] 



■ y^cont iT^L' 

sou,+ sys,± 

cont _L ^[cal] ^^ 



usually unable to determine Tsou,* ■ ^ sys,± 
— required for Eq. (EU and dS} — but only T^^^^^ -h r,y, + 
a consequence the inferred calibration is wrong at frequencies 
where the spectral-line ghosts appear. This does not only affect 
the amplitude of the ghosts, but even their shape gets distorted. 
To show the effect of this, we use a simple model in which differ- 
ent source intensities (Jsou = 5, 10, 20, 50, 100 K) were added to 
a hypothetical system temperature, Tsys = 100 K (top panel). 
The bandpass was assumed to be flat. For visualisation pur- 
poses, the spectral lines are triangular-shaped. Increasing the in- 
put source temperature with respect to Tsys decreases the line 
ratios of Tghost to Tsou after applying frequency shifting (mid- 
dle panel). In addition the distortion becomes more and more 
pronounced. Finally, in the bottom panel the result of a simple 
shift-flip- sign-averaging (so-called folding; see also Eq. ([37])) 
is presented. The legend contains the reconstructed peak tem- 
peratures, which are systematically deficient for larger ratios of 
Tsou/Tsys- This folding algorithm should clearly be avoided for 
cases where Tsou ^ T^y^. 



there is a residual additivd^ baseline depending on the relative 
changes in Tgys and Tcai with frequency. These can make the fi- 
nal baseline fitting a challenging task, especially in the presence 
of standing waves (SW). 

One thing to note is that usually we are unable to determine 
+ T^^^^^ -\- Tlyfl and can only compute 



+ r. 



[cal] _ 7-line 



sys,+ 



_ T^ll 



^ This is an important point, as it allows for a simple subtraction of a 
baseline fit, without destroying the calibration. 



^s^ou,± + ^sys!± which then is substituted into the right-hand side of 
Eq. (l34l) and (1351) . This is an essential drawback. Consequently, 
the estimated 7'sig,ref are wrong for all frequencies where T^JJIf (y± 
Ay) ^ Tsys is not fulfilled (see Fig.fT4l). 

Luckily, this usually only happens for frequencies at which 
the spectra l-lin e ghosts appear (except for complicated cases; 
see Section l46l) . If one uses both Eq. ([34b and ([35]) and applies 
shift-and- averaging as described above, this is not a problem. 
However, if just one of these two equations is used and instead 
a so-called fold' procedure (e.g., the Class task fold from the 
GiLDAS packago is applied, one will achieve incorrect results if 
Tsou ^ ^sys (Fig-IHl bottom panel). This is not a rare case, espe- 
cially for observations in the decimetre- wavelength regime. The 
term folding denotes the use of spectral-line ghosts themselves 
to add to the positive emission line, by changing the sign of and 
both shifting and averaging the signals 

Uf^^\y^Ay)-f^^\y-Ay)]. (37) 



4.2. Setting up simulations 

The input spectra for frequency switching are similar to those 
shown in Fig.[Tl except that each measurement phase contains 
the Tsou contribution. They are also shifted according to an LO 
frequency of iAvLo = ±10 MHz before multiplying by the 
bandpass Gif (Grf = 1). A different standing wave frequency, 
Vsw = 10 MHz, is also used for the method that we describe in 
Section [44l (since vlo = 10 MHz must be a multiple of Vsw for 
this second method). Fig. [15] shows the resulting input spectra as 
they would have been measured by the backend. We note that for 
frequency switching the measured sig and ref phases belong to 
different RF frequencies. To clarify whether an operation is per- 
formed in either the IF or RF domain, we assume a hypothetical 
IF centre frequency of 150 MHz and plot figures accordingly. 



4.3. Using the two switching phases individually to obtain 

^ sys 

As in the previous sections, the signal of the calibration diode is 
used. Computing 
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(38) 



(39) 



and inserting these expressions into Eq. ([34] ) and Eq. ([35] ) leads 
to 



y^[eal] _ 
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ref 



(^Hne,+ + ^ref 7'cal,+ + [^eal,+]) ' 
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(40) 



(41) 



One needs to neglect the line emission contribution, Tsou, on the 
right hand side of the equations, but this has no serious impact, 
except that spectral-line ghosts have unexpected amplitudes. 

Fig.[T6]displays the a:^.^ ^^^ spectra, along with appropriate fit- 
ting models. Using these models leads to the calibrated spectra 
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Fig. 15. Shifting the input spectral densities with iAvLo = 
±10 MHz and multiplying by the (IF) bandpass shape (see 
Fig.O returns the 'measured' spectra /^sig,ref- The upper panel 
shows the results for the simpler case, while the lower panel 
includes a standing wave contribution (vsw = 50 MHz). Note 
that for the method described in Section l44l we used a smaller 
standing wave frequency (Vsw = vlo = 10 MHz) in contrast 
to the other examples in order to avoid resonances (see also 
Section l4AT]) . 



shown in Fig. [T7] (top panels). Whilst the more simple case with- 
out SW shows no residual baseline owing to the specific form 
of Eq. ([36b in the resulting spectrum, the SW case displays a 
strong residual pattern. For the lower right panel in Fig.[l7j we 
subtracted a sine- wave signal modulated with a polynomial. 

4.4. Using both switching phases together to obtain T^ys 
We parametrise 
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(43) 



using only the continuum contributions to compute a model (see 
Fig.[T8]). The SW frequency was set to 10 MHz, for reasons dis- 
cussed in Section l44.1l 



The SW case requires a relatively complex model of the form 

m(y) = [L(n, v - Ay) - L(n, v -\- Ay)-\- 

A sm(a(y - Ay) -\- b) - A sin(a(y -h Ay) -h b)] x (44) 

[P(n, y -\- Ay) -\- A sin(a(y -\- Ay) + Z?)]"^ , 

where L(n, v) is a polynomial function of degree n. Fig.fTSlshows 
the results of these fits. In the presence of SW (right panel), one 
must be sure to have a good initial estimate of ygw in order to 
obtain a meaningful result. This may impose difficulties in prac- 
tice, though it is likely that using a particularly long integration 
time (or using several data sets) will aid the measurement of v^^. 
From Eq. (l42l) and Eq. (1431) , we obtain 
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(46) 



The term Arcai,±/(/ - /^^^) contributes significantly to a andyS 
and must not be neglected. This correction term is the only fun- 
damental diff'erence between the frequency and position switch- 
ing equations, which otherwise are very similar. This is because 
the ref phase is not independent of the sig phase, but contains 
only frequency- shifted quantities, introducing a degeneracy. As 
for position switching, a problematic situation occurs if / ~ /^^^ 
since the denominator in Eq. (1451) and (l46l) becomes singular. 
While in the case of SW this is a frequent phenomenon (see 
Section l4AT]) , it may also occur when Ar^'^^J^^ ^ 0, Arsys,± ~ 0, 
and Arcai,± ~ 0. At first glance, this might appear to be rather 
unlikely, but it can indeed arise if the continuum and Tcai slopes 
are not very steep in relation to the LO shift. 

Using the substitution v ^ v- 2 Ay in Eq. (l45l) and (l46b , one 
can also infer 
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(47) 
.,y,,_=Av-2Ay). (48) 

Subtracting Eq. (gS]) from (|471) and Eq. (gS]) from diU results in 
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p{v-2^y)-|3{y) = rZ. 
Inserting Eq. (|49|l and dUll into Eq. (|33ll leads to 
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(42) Similarly, 
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Fig. 16. The spectrum of k^^^ ^^^ for the simpler (left panel) and SW (right panel) cases and the inferred fitting models. Spectral 
windows are set around parts containing the spectral lines (not clearly visible in the plots due to noise), otherwise the models would 
have been influenced to some extent. 
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Fig. 17. Using Eq. ([36l) and the fitting models from Fig. [161 one obtains the calibrated spectra. The upper left panel shows the result 
for the non-SW case, and the upper right panel for the SW case. A baseline was subtracted for both in the lower panels. In the 
SW case, the residual baseline is very complicated. The continuum of the source is not reconstructed. 



Finally, to average all four phases to form the final spectrum, using Eq. (l45l) , even if rcai(v) is known. This is due to the degen- 
each resulting spectrum must be shifted to achieve that the emis- eracy described above. 

sion lines lie at the same position in the spectrum. One advantage of this method, despite being more complex, 

is that on the left hand side of Eq. ([52b to Eq. ([55]) , there is 

We note that, in contrast to position switching, one cannot no continuum contribution (in contrast to the method presented 
reconstruct T^^^^ and Tgys as separate quantities, only T^^^^ + Tgys in Section [43]) . This means that SW do not propagate into the 
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Fig. 18. To correctly reconstruct the original fluxes, our proposed method utilises models of the intermediate spectra (P^^f^ 



^'ef )/^'ef (see Eg. (l42l) and (143))). Appropriate windows around the spectral lines must be set as one is interested in the con- 
tinuum contribution only. The left panel shows the results for the simpler case using second-order polynomials that are able to 
describe the baseline. For the standing-wave case (right panel), a more complicated fitting model m(y) had to be applied (see 
Eq. dH). 



final spectrum producing complicated baselines as in Fig.fTTl 
However, SW can still have a negative influence on the model 
fitting such that the calibration could be aflfected. This means 
that if one has to reduce frequency switched spectra having a 
SW contribution, this method may be the most suited in terms of 
residual baselines, however the calibration accuracy is lowered. 
Results for the non-SW and SW cases are displayed in 
Fig. [19] (upper left and upper right panels respectively). The 
results are less accurate than for the method presented in 
Section |43l Nonetheless, further investigation, e.g. studying the 
distribution of the flux calibration error using several hundreds 
of simulated spectra (see Section [6^ shows that the calibration 
is still unbiased. The error distribution is however wider than for 
the other methods, owing to the use of two independent fitting 
models and the relative errors between both fits then having a 
rather large impact on the final results. The method is even less 
robust for the SW case. 

4.4.1 . Problems caused by standing wave contributions 

During our simulations, an important issue emerged. Inspection 
of Eq. (|45]) reveals that a(y) can become singular, if / = /^^^ or 



^ cal,+ 
Tcal- 



^ cont,+ "•" ^ sys,- 



ysou 
cont,- 



+ ^ sys,- 



(56) 



Under normal circumstances, this condition will not be ful- 
filled. However, in the presence of a standing wave contribution. 



4.5. Determination of T^ys neglecting frequency 
dependencies — the 'classical' approach 

On the basis of Eq. (1381) and Eq. (|39l) , but explicitly not restrict- 
ing ourselves to continuum contributions, we find that 



^ line,- "^ ^ cont,- "^ ^ sys,- 
^ line,+ "^ ^ cont,+ "^ ^ sys,- 



f cal,- 






(P',:l + Pref)-(P'r^}-P.ef) 



SlgJ 

cal 
ref 



^(Pr^-Pr.f) 



(57) 
(58) 



Neglecting line-emission contributions completely and treating 
Tsys and Tcai as constants, one can write 



^ cont,- ^ ^ sys,- - \^ cont,- ^ ^ sys,-/y 
Dcal 






sig/y 



cont,+ "•" ^ sys,+ - \^cont,+ "•" ^sys,+ /y 



= T, 



cal 



(P'3^PM)-{Pff-PM\ 



2{P'S-P,,f), 



ref 



(59) 



(60) 



Under the assumption of frequency independence, Eq. (|34] ) and 
Eq. (|35]) simplify to 



rsw,+/- to rsys,+/-, there will be several frequencies in the spec- Tgou,- - 7^sou,+ = (7'Hne,+ + ^cont,+ + ^sys,+ [+7^cai]) 



trum where Eq. (1561 ) is true. This results in 'resonance' spikes in 
the reduced spectra, shifted according to vlq) (see Fig.[20l upper 
right panel). These occur if the LO frequency is not (close to) a 
multiple of the standing wave frequency. The middle panels of 
Fig.[^show (P^f ^^ - /"l^f )//'|ef ^^^ the case where the LO shift 
matches Vsw (left) and does not match Vsw (right). The lower pan- 
els depict Eq. (l56b , i.e. the resonance condition in terms of tem- 
peratures. Consequently, to avoid this problem one must ensure 
that rsw,+ ~ T'sw,- by matching the LO shift to the frequency of 
the standing wave, i.e. vlq ~ ^Vsw, ^ e IN+. 



-^ sou,+ -'sou,- - (Mine,- "*" -'cont,- "•" -'sys,- L+-' calJj 



p[cal] _ p[cal] 
sig ref 

p[cal] ' 

ref 

(61) 

p[cal] _ p[cal] 
ref sig 



p[cal] 
sig 



(62) 



Figure [2n contains the resulting spectra after applying the 'clas- 
sical' method, shifting, averaging, and subtracting a baseline. 
The recovered spectral line intensities are systematically wrong. 
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Fig. 19. Inserting the models f^^^^\y) into Eq. ([52b results in a unbiased flux calibration, as shown in Section[6]2j The figure shows 
the outcome for the non-SW (left panels) and SW (right panels) scenarios with (lower row) and without (upper row) baseline 
subtrac tion. It appears that the method is quite significantly aflfected by small errors in the fitting models of / and /^^^ (compare with 
Section [6^ . Furthermore, to obtain acceptable fits, we had to choose very good prior model parameters, e.g. from a low-noise run. 
Nonetheless, compared to the method in Section |43l the SW do not directly convert into a sine- wave imprint in the final spectrum 
(compare upper right panel with Fig.[T7]l, although relatively often higher order eflTects are caused when the model fits are inaccurate. 



In the SW case, the baseline has a complicated pattern, which 
can hardly be described by a simple model. 

We note that simple spectral averaging to calculate Tsys is 
aflTected by spectral lines and spectral-line ghosts, introducing 
even more inaccuracies. 



4.6. The realistic case — limits of frequency switching 



In Section [X6l we have shown that position switching also works 
well for more complicated examples, where the system temper- 
ature and calibration signal are fluctuating. The question now is 
whether one can also apply the frequency switching schemes to 
this present example. 

For the second calibration model, it is usually impossible 
for the reasons discussed in Section l44.1l Fig. [22] illustrates 



T,,i,JT,,i,-^ndmZ^^T.^ 



- sys,- 



-)/(-^cont,- "•" -* sys,- 



). Each intersec- 



tion of the two functions would lead to a singularity in the final 
spectrum rendering a fair fraction of the spectral range unusable. 
Furthermore, as we have seen in Section l44l the accuracy of the 
second method is generally lower than that of the first. 

Does the first method (see Section |43]) work? In short, the 
answer is yes, but the result is still unusable (see Fig. [23]). While 
the reconstructed spectrum is correctly calibrated, the residual 



baseline (which is also correct) has such a complex shape that 
one can hardly distinguish the spectral line from the structures in 
the baseline. This has nothing to do with our calibration scheme 
but is an intrinsic problem with frequency switching. Inspecting 
Eq. ([36b reveals that not only does the source produce spectral- 
line ghosts, but there is the term 

Tl^f\y) - \Tl^f\y + 2Ay) - ^Tl^lHv - 2Ay) . (63) 

Hence, as soon as the system temperature (and Tcai) has spectral 
features on the same scale as the astronomical lines, frequency 
switching will fail (or rather the astronomer who has to extract 
useful information out of the spectrum). This is an important 
aspect to be considered for future broadband instruments. 



4. 7. Combining position and frequency switching to handle 
^RF,- ^ ^RF,+ cases 

At the beginning of SectionlH we discussed how frequency 
switching is limited to cases where Grf- = Grf,+ . This is un- 
fortunate because it substantially limits the number of possible 
applications as almost no real receiving system will provide a 
constant RF gain, at least if moderate or large bandwidths are 
involved. 
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Fig. 20. If the LO shift vlo is not a multiple of the standing wave frequency Vsw, the quantity / - /^^^ can have several roots causing 
'resonances' in the final spectrum. The left panels have Vsw = 10 MHz matching vlo = 10 MHz, while for the right panels we chose 
Vsw = 50 MHz. The upper panels show the resulting spectra (simulations were performed without a noise contribution for a clearer 
visualisation). In the upper right panel, spikes do appear. The middle panels show the model fitting through (P^^^^^ - P^^^f^) / P^^^f\ 

While in the left middle panel /^^^ > / for all frequencies, in the right middle panel several crossings occur. The lower panels 
depict Eq. (l56l) , which represents the condition for resonances to appear. 



One may however combine a frequency- and position- 
switching measurement, preferably of a calibration source, to 
infer the term Grf,-/Grf,+ . This is sufficient to open the fre- 
quency switching method to a broader range of applications 
where Grf,- 9^ Grf,+, as we show below. Without any restric- 



tion with respect to Grf, Eq. (1331) reads 



3[cal] 



^[cal] 
ref 
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+ r. 
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Fig. 22. Applying the second calibration method is impossi- 
ble for more complicated system temperatures and/or calibra- 
tion signals. Using the values from Fig.[TOl the two functions 

rcai,+/rcai,- and (TZt,^ + T,y,,^)/(TZt,- + ^ys,-) have several 
intersections, which would lead to singularities in the result. 
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frequency switching and the first calibration method (see Fig.fTUl 
for temperature inputs). The only diff'erence in comparison to 
the position switching example is the use of a single Gaussian 
with an amplitude of 10 K. Despite the larger line amplitude, the 
emission peak is not very pronounced in the result. 
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where the term T^o^^+fT^o^- is known if a calibration source is 
used. A simple rescaling 
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ref 



frequency- switching equations is sufficient to fully incorporate 
the Grf,- 9^ Grf,+ cases. 

The On-Off measurement of the calibration source has of 
course to be repeated on appropriate timescales to account for 
time-dependent instabilities in the RF gain. However, in most 
cases this will allow a much longer effective integration time on 
source than a position-switching measurement, which in the past 
has been considered the only viable option for 'futile' cases. 

4.8. Least squares frequency switching 

Least squares frequency switching (LSFS) is a technique devel- 
oped by Heiles (2007), which makes use of not only two, but sev- 
eral different LO shiftfl. The underlying mathematical problem 
can be approximated by a linear equation and efficiently solved 
using maximum likelihood methods. A great property of LSFS is 
that it does not produce spectral-line ghosts. Hence, this method 
is the only option among frequency switching techniques for 
cases where blending of emission lines with spectral-line ghosts 
would be a problem, or where the system temperature or Tcai 
has structure (see also Section l46] ). Improvements to the orig- 
inal LSFS method were proposed by Winkel & Kerp (2007) to 
incorporate RFI treatment and to increase its robustness in the 
presence of strong emission lines. 

In contrast to classic frequency switching, one is able 
to determine the IF bandpass curve, as well as the product 
Grf (^sou + ^sys)- The latter however shows that as for normal 
frequency switching, only the IF gain can be dealt with. Hence, 
if Grf is frequency dependent (which is often the case especially 
for large bandwidths), frequency switching will lead to calibra- 
tion errors. A similar method to that discussed in Section I^TTl mav 
be applied to work around this problem, i.e. combining LSFS 
with a position switch. It is then possible to even disentangle 
Grf and Gif, which can be considered a great advantage over 
classical frequency switching. Some overhead would clearly be 
involved to conduct (repeated) necessary On-Off measurements 
toward calibration sources. 



5. Direct determination of the bandpass curve 

For the sake of completeness, we now present a third method 
that may be suitable for obtaining a sufficient calibration. It uses 
averaging over a large number of sequential spectral dumps, P/, 
in order to measure the bandpass directly 



t^i - fi - CtrfCtif [I- + 1 i^yg - 1^ -I i,sys j 
= GrfGif [T^cal + AJsou + AJsysJ . 



(71) 



The correct rcai(v) clearly needs to be known here. Since AJsou + 
AJsys ^ Teal, as continuum levels between two subsequent 



^ This might require changes to the hardware of a receiving system to 
allow for more than the usual two LO frequencies. We note that LSFS 
does not deliver calibrated data but only normalised spectra that must 
be flux-calibrated using other measures (e.g. using Tcai) later on. 



dumps within one scan can be assumed to be approximately 
equal, then 
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sysl , 



(p-'-P>^ 



^cal 



= r(v). (72) 



The function y(y) is of course subject to noise, which owing to 
averaging is however much smaller than in the raw data of a sin- 
gle dump. One may even apply spectral averaging (smoothing) 
to further decrease the noise in y(y). It follows that 
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cal- 



(73) 



An important difference from the position switching method is 
that one can only infer Tsou + ^sys, hence afterwards one can- 
not distinguish the background and source continua. The back- 
ground continuum also varies with time, and thus the baseline 
fitting must be carried out on appropriate timescales. A second 
drawback is that the bandpass/gain must not vary over the course 
of the observation, otherwise the procedure would only return an 
average bandpass such that the individual spectra are wrongly 
calibrated. Whether this method is applicable depends on the 
properties of the receiver. Nevertheless, the performance of the 
direct bandpass estimation may still be better than the use of fre- 
quency switching (one also has full-time on-source, in contrast 
to position switching). 

This method is us ed for the Effelsberg-Bonn Hi Survey 
(EBHIS; IWinkel et al.l 120 IQ) , where neither position nor fre- 
quency switching is applicable. For the 21 -cm line at galactic 
radial velocities, there is simply no suitable Off position as there 
is emission in every direction. The receiver also exhibits a com- 
plex multi-modal SW pattern, which would be hard to model 
using analytic fitting functions. Frequency switching is impos- 
sible as the receiver has a signifi cant bandpass ripple in the RF 
part. For more details, we refer to IWinkel et al.l (|2Q1Q[ l2Qlll) . 

A simulation was conducted using the same set-up as for the 
position switching, with the only difference being that no Off 
position was needed. In Fig.[24l (upper left panel), the input spec- 
tra after multiplying by the bandpass curve are shown (including 
SW). 

Figure|24] also presents P^^^ - P for a single spectral dump 
and (P^^^ - P}t, as computed using the average of 1000 spec- 
tral dumps, the latter showing significantly less noise. The upper 
right panel presents the reconstructed spectrum. In contrast to 
position switching, the full continuum contribution Tsou + Tsys is 
obtained (including SW). Therefore, to extract the pure spectral- 
line contribution one has to compute an appropriate baseline (in 
our case the superposition of a polynomial and a sine wave). 
The resulting spectra (after baseline subtraction) are shown in 
the lower panels. We use Tca\/Tsys,oS ~ 15% (left panel) and 
Tca\/Tsys,ofi ^2.5% (right panel). In the resulting spectrum, the 
noise is of course higher in the latter case. It is however much 
harder to obtain a noise-free description of the bandpass curve 
(by filtering or model fitting) than in the former case. The ro- 
bustness of the method therefore also depends on the ratio of 

^ cal to 1 sys,off • 



6. Flux density calibration error distributions 

We analyse the spread in the recovered flux values using the 
three diflferent calibration methods. The distribution of flux er- 
rors is a good indicator of the robustness of the calibration 
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Fig. 24. Top left panel: Input spectra (including SW) used for the direct method to determine the bandpass curve. It also contains 
peal _ p for a single spectral dump and {P^^^ - P)t as computed using the average of 1000 spectral dumps, both scaled by a factor of 
five for improved visualisation. The upper right panel shows the reduced data after applying Eq. (1731 ). Note that the standing wave 
contribution is included in the result. The lower panels show the final spectra after baseline subtraction for Tc^i/T^y^^of^ ^15% (left 
panel) and rcai/rsys,off ~ 2.5% (right panel). 



schemes and also reveals potential biases. The exact values ex- 
pected depend strongly on the specific choice of the various tem- 
perature contributions (especially the slopes), bandwidth, and so 
on. Our analysis should therefore be treated as a qualitative study 
and not a quantitative analysis. 

To estimate the distribution of flux errors, we used the sim- 
ulations (of the more simple example presented in Section [3Jl) 
as described in the previous sections, but repeated each of them 
one thousand times, each time with a difl'erent noise realisation 
but all other parameters constant. In this way, we were able to 
explore the stability of the various methods. 



for frequency-dependent system temperatures results in an un- 
biased estimate of the true flux values. Interestingly, the classi- 
cal method has a similar scatter (cr) in the derived flux values. 
Nonetheless, intensities are systematically wrong in an unpre- 
dictable way. In the absence of SW, the two more sophisticated 
methods provide similar results. This argues in favour of the first 
method relying only on one fitting model, while the second re- 
quires two. In the presence of SW, the second method is still un- 
biased but has a somewhat larger scatter. We conclude that the 
calibration scheme described in Section [T2l is optimal for posi- 
tion switching measurements. 



6.1. Position switching 

The resulting histograms of the relative flux error for the dif- 
ferent cases are shown in Fig. [25] for each of the three simu- 
lated Gaussian emission lines. The left (right) panels show the 
case without (with) standing waves (SW). The upper row con- 
tains the results provided by the method described in Section [T2l 
(using only the Off position), the middle panels those for when 
both the On and Off positions were used (see Section [33l) , and 
in the lower panels the classical calibration scheme was applied 
(Section [T4l) . The bin width was the same for each histogram, 
to enable the amplitudes (normalised number counts) to be di- 
rectly compared. We fitted a Gaussian to each distribution of er- 
rors, the parameters (mean value yu, width cr) of which are stated 
in the plots. Using the proposed calibration methods to account 



6.2. Frequency switcliing 

As above. Fig. [26] illustrates the results for the three diff'erent 
calibration methods for frequency switching. The proposed cali- 
bration schemes provide unbiased intensity estimates, while the 
classical method is aff'ected by systematic errors. Unfortunately, 
the two-model approach (middle panels) results in a much higher 
spread around the true value than for the position switching sim- 
ulations, as the diff'erence between / and /^^^ is smaller — a di- 
rect consequence of frequency switching. As for position switch- 
ing, the first method (see Section l43l) is generally more suitable, 
as the associated error distribution is narrower. Nevertheless, the 
second method, Section l44l might be considered if baselines be- 
come a problem, as the latter approach is less aff'ected by the 
standing wave shape in the final spectra. 
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Fig. 25. To estimate the robustness of the different methods, we repeated each simulation a thousand times. The resulting flux 
errors for the three Gaussian components were evaluated. The histograms show the error distribution for the two cases, left (right) 
panels without (with) SW. The upper rows show the results for the method described in Section lT2l using only Off positions to 
estimate Tcai/T^ys- The middle panels were produced using the procedure from Section l33l The lower row contains the result for 
the classical method (T^y^ = const, and Tcai = const.). The bin width has been kept constant for each histogram, such that the 
amplitudes can be directly compared. For each distribution of errors, we fitted a Gaussian distribution, the parameters (mean value 
ji and width cr) of which are stated in the plots. 



We note that we use the same noise contribution in both the 
position and frequency switching methods. Owing to the in-band 
frequency switching, the resulting noise in the final spectrum is 
smaller by a factor of ^/2. Whether the lower level of noise can 
compensate for the potentially lower quality of the frequency 
switching results is a decision that has to be made by the ob- 
server, because it depends on the chosen instrumental set-up (e.g. 
receiver, SW) and observing scheme. 



7, Measuring Tcai 

The reliability of the calibration in all the methods that we have 
described strongly relies on accurate knowledge of the noise 
diode's spectral density. Consequently, some eff'ort should be in- 
vested in the measurement of Tcai. 



7.1. Using hot-cold calibration 

The only means of determining Tcai directly is to apply a hot- 
cold calibration. This is achieved by pointing the feed horn to an 
absorber at ambient temperature and subsequently to a tank con- 
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Fig. 26. As in Fig.[25l the plots show the error distribution for a thousand simulations after applying frequency switching. For the 
first method, Section l43l the resulting histograms are somewhat smaller, which can be attributed to the slightly smaller residual 
baseline RMS (~ 80 mK for frequency switching, ~ 120 mK for position switching) in our simulations. The situation is much worse 
when using the second calibration scheme, Section l44l The distributions are broader, even more so in the SW case, owing to the 
relatively small diff'erence between / and /^^^. 



taining liquid nitrogen at 77 K. The experiment is time critical as 
after a while the feed itself starts to physically react to the low 
temperatures causing instabilities in the system. 

Since the temperatures of the hot and cold phases are known, 
one can easily compute the transfer function (i.e. bandpass curve 
and calibration factor) that converts the input spectra (in K) to 
spectrometer output (in arbitrary units) 



G(y) 



r,[ca\] 
hot 



r,[cal] 
cold 



T--[cal] 
^hot 



T-'[cal] ' 
^ cold 



(74) 



using Thot = 300 K and Tcoid = 77 K. Dividing the input spectra 
P by G, one obtains calibrated data. Finally, after calibrating the 
spectra we can compute the contribution of the noise diode Tcai 



by subtracting the non-cal from the cal phase 



^ cal - ^ hotxold ^ hot,cold 



peal 
hot,cold 



hot,cold 



(75) 



In Appendix lB.il we provide an example of a hot-cold calibra- 
tion measurement of the 6. 5 -mm receiver at the 100-m telescope 
Eff'elsberg. 

Since a small error in the determination of Tcai has a large im- 
pact on the later data reduction, hot-cold measurements should 
be performed with great care. There are some potential issues, 
which we briefly discuss in the following. 

Firstly, the hot load at room temperature, i.e. ~ 300 K emits a 
much larger amount of power than the sources usually observed. 
This diflTerence of 20 dB should not be a problem for any modern 
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receiver and backencfl. However, it is important to level the sys- 
tem accordingly, otherwise the amplifiers might be over-driven 
on the hot load leading to compression of the signal. 

Feed horns in the primary focus of a telescope have a rather 
large beam width because they need to illuminate the antenna 
dish. As a consequence, they have large spill-over angles, which 
require the hot and cold load to be large enough such that all side 
lobes of the feed also point towards the load. To achieve this, one 
should place the load, usually an absorber, in a metallic tank and 
plunge the feed into this as much as possible. For the hot load, 
side lobes play a lesser role since the environment is usually at 
the same (ambient) temperature as the absorber, e.g. when the re- 
ceiver is in the lab or in the focus cabin of the 100-m telescope. 
However, if the receiver is not in a closed environment, radia- 
tion from the sky may be reflected and enter the horn, changing 
the net radiative temperature of the load. Likewise for the cold 
load, radiation from the environment must be avoided to as great 
an extent as possible. In practice, the radiation temperature of 
the cold load is up to 2 K or 3 K higher than the n ominal 77 K 
(whic h itself depends on the atmospheric pressure, lUlich et aP 
119801) owing to this leakage. 

The cold load is tricky to handle. The absorber must be com- 
pletely immersed in the liquid nitrogen, the amount of which 
can be a logistical issue at a remote site or for large receivers, 
e.g. multi-horn systems. Furthermore, both the vaporisation of 
nitrogen, and frost that can potentially form on the feed enforces 
strict time constraints. Frost reflects part of the incident radiation 
and changes the properties of the feed. If several hot-cold mea- 
surements are desired for one receiver, e.g. when the receiver 
has several sub-bands, one should switch between hot and cold 
regularly while ensuring a constant height above the absorber. 

The absorption material itself must fit the observed fre- 
quency. For low frequencies, below several 100 MHz, it can be 
a problem to even find an efl&cient absorber, which should also 
have a good heat conductivity in order to reach the temperature 
of liquid nitrogen suflficiently quickly. 

Unfortunately, calibration diodes are known to drift with 
time. It is therefore essential to repeat hot-cold measurements 
at least every few months if high precision is desired. 

7.2. Using continuum sources with well-known flux 

Another way of obtaining a good estimate of the Tcai spectrum 
is to utilise strong continuum sources with a well-known flux 
density. For example, using position switching one can easily 
invert Eq. ([T6l) 



Teal = 



Kof^- 



-Ppff 



PqS 



(Koff + 1) 



r)cal_ 



■^off 



(76) 



Depending on the S/N of the recorded spectra, it may not even 
be necessary to use a model of Kq^. This of course increases the 
noise in the Tcai measurement somewhat, but has the advantage 
that the result is independent of the choice of a specific model. 
We present an example measurement using the 6-cm secondary- 
focus receiver at the 100-m in Appendix lB.2l 

The true brightness temperature, r*^^ = ^ap ^sou exp (tqAM), 
(or rather the associated flux density S = Ik^T^^JA, where A 
is the physical area of the dish) is usually tabulated in a contin- 
uum flux catalogue, which is of course not the measured antenna 



^ At least for the FPGA-based backends that are in use at most major 
facilities today. Older devices, such as auto-correlators with their low 
input quantisation, might cause problems. 



temperature Tsou- As a consequence, one not only needs to take 
into account atmospheric dampening but also the aperture ef- 
ficiency, 7/ap. To determine z/ap for a moderately- sized or large 
single dish telescope, one has no choice but to measure it using 
astronomical (point-like) sources. This is usually done using a 
hot-cold measurement (to absolutely calibrate Tcai, which yields 
Tsou after application of Eq. (l76l) ) and — in a timely manner — 
relate it to the known flux density of a calibration source. Thus, 
how is the latter determined in the first place? The brightest ra- 
dio sources in the sky have been observed with well-calibrated 
horn antennae (or a dipole arrays), the response of which can 
be theoretically calculated or precisely measured (since they are 
smaller). This served to establish the (frequency-dependent) flux 
values for sources such as Cas A, VirA, or CygA. From these 
primary calibrators, flux scales for additional sources were de- 
rived using relative measur ements. For a detailed discussion , we 
refer to the literature(e.g., lFindlaylll966i iBaars et aDll977l and 
references there-in). It is clear that if the use of continuum flux 
catalogues is desired, the r/ap values should be regularly updated, 
using hot-cold calibration. 

Furthermore, to use a continuum source for a measurement 
of Teal the condition AJsys = must be fulfilled (see Eq. ([8])). In 
practice, this can be an issue, for example the SW contribution 
can depend on the incident continuum radiation (I Winkel et al.l 
2011). The observer should at least ensure that the On and Off 
measurements track the same elevation angle(s), such that the 
amount of spill-over (ground radiation) and the atmospheric con- 
tribution are similar. 



8. Summary 

When dealing with position- or frequency- switched observa- 
tions, we have shown that in many cases it is of utmost im- 
portance to account for the frequency dependence of the system 
temperature and noise calibration signal in the flux calibration 
task. Simplification of the problem, i.e., treating Tsys and Tcai as 
constants, leads to a systematic calibration bias, which can have 
a strong impact on the scientific results from reduced data. We 
have proposed several methods to circumvent this bias at the cost 
of a greater complexity in the data reduction. 

The conclusions that we can draw from this paper include: 

1 . In general, position switching provides much more accurate 
results than frequency switching because it is less aflfected by 
a variety of potential issues, such as standing waves. If the 
system temperature and/or calibration signal have strongly 
fluctuating spectral features on similar or smaller scales than 
the astronomical signal that is to be observed, frequency 
switching is not applicable at all, because of the spectral-line 
ghost issue. In these cases, least squares frequency switching 
might well be worth the higher technical and computational 
effort. Unfortunately, 50% of the observing time is efl'ec- 
tively lost when using position switching, such that the final 
choice between the two techniques is reduced to the question 
of whether a longer integration time on- source compensates 
for the potential problems. The answer to this question is 
of course diff'erent for every receiving system and observ- 
ing condition. For non-constant Grf(v), one cannot achieve 
a correct calibration using frequency switching, but we have 
shown how a combination of position and frequency switch- 
ing can be applied in these cases. 

2. Another significant advantage of position switching is the 
ability to reconstruct the continuum flux of the source. This 
also makes it possible to use known continuum calibrators 



21 



B. Winkel et al.: Unbiased flux calibration methods for spectral-line radio observations 



to determine the Tcai spectrum, which is essential to track 
time-dependent drifts. 

3. The presence of standing waves make sophisticated calibra- 
tion more difficult, as the quality of a non-linear least squares 
fitting involving sine-functions relies heavily on the initial 
model parameters. This problem is even more pronounced 
for the second approach, which makes use of two indepen- 
dent functions (see Sections [33l and l44l) . The use of signal 
filters, instead of parametric models, to suppress noise, e.g. 
in a:" \ is a solution to this problem but at the cost of intro- 
ducing an additional (correlated) small noise component. 

4. The Teal intensity should be sufficiently high. All proposed 
methods utilise the diff'erence (and its inverse) between the 
spectra that contain the noise diode's signal and those that do 
not. A higher S/N helps us to avoid numerical instabilities. If 
Teal is too large, the average system temperature is of course 
much higher than necessary. Usually Tcai/T^ys ~ 20% is a 
suitable value. 

5. Although under certain circumstances, the sophisticated 
methods presented are not as robust as the classical method, 
e.g. in the presence of standing waves, they allow for an un- 
biased calibration that is of utmost importance to all scien- 
tific investigations. 
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perature can be modelled (e.g., using one of our methods) and 
hence is noise-free, we can rewrite Eq. (|9]l 



^ sou + ^^ sys — (^ sys Qffjmodel 



/y.[cal] X 



sys, 



'model 



p[cal] _ p[cal] 

p[cal] 
^off 

y.[cal] 
-^ on _ 

^[cal] 
^off 



(A.1) 



where we have used the shorter notation T^^^^ = tII^^ -\- Tlyf^^ 



and T 



[cal] 



= T 



[cal] 



, ^^ _ ^^ ^^Q. The latter two quantities are affected by 
noise according to the radiometer equation (see Eq. (1111). The 
resulting RMS is then 



^[cal] _ /7-[cal] 



_ /T^LCciij X 
rms — y^ sys,off^niodel ' 



/y.[cal] x2 
v^ rms, on/ 



/y.[cal] x2/r 
^^rms,off^ ^^( 
/^[cal]x4 
^^off ^ 



[cal]x2 



A 



/y.[cal] 
v^ rms, on 



^2 ■ /^[cal] \2 
^ ^ y^ rms,off^ 



/^[cal]x2' 
^^off ^ 



(A.2) 



which is at least a factor of ^/2 larger than the noise in each 
individual phase (if the continuum of the source is negligible). 
Likewise, the S/N is a factor of ^/2 smaller because in position 
switching two raw spectra contribute to the noise power but only 
one of them contains a signal. This is the trade-off one has to pay 
when applying a switching technique. The final noise level after 
averaging both phases, Taver = (^sou + T^f^ll, is then 



-\4 



(r™s)2 + {T\ 



■' sys,off 



»r^ 



V2' 



(A.3) 



which is roughly a decrease by ^|2 relative to each of the 
cal/non-cal phases individually. This implies that the final RMS 
is only about the same level as the (mean) noise in the individual 
phases! As discussed in Section [T2l one possible way of improv- 
ing on the final noise would be to smooth the reference spectrum 
prior to division, when only the noise of the on phase would then 
propagate into the final result. 

There might be cases, where rsys,off ^ T^^^ and Eq. (IA.3b is 
dominated by T^^ increasing the noise in the averaged spectrum. 
The following weighting scheme could then be used to improve 
the situation 



,cal T-'cal 



wTsou + ^^l 



(A.4) 



where (h'['=^")-i oc var[r['=''"] = (rJnfsV. One can easily infer that 



real 



(A.5) 



Appendix A: Theoretical noise levels for position 
and frequency switching 

A.1. Position switching 

Using Eq. Q, we can calculate the theoretical noise level, Trms, 
of the cal and non-cal phases. Assuming that the system tem- 



which leads to smaller noise values in the averaged spectrum 
than Eq. (IA.3I) . In the limit rsys,off » T^^^, both equations return 
the same result. 

A practical problem, however, is to estimate the noise spec- 



r[cal] 



trum r^s (^) needed to determine the weighting spectra w^^^\v) 
with sufficient precision, which can be complicated in the pres- 
ence of either emission lines or radio frequency interference. 
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yA.2. Frequency switching 

As in the case of position switching, we modify Eq. (l34l) and 
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-1 



(A.7) 



where T^f = Tsou,- + Tl^fl. and 7^^" = Tsou,. + tI;1% We 
note, that for the modelled (effective) system temperature we 
assume that only continuum contributions could be inferred, as 
discussed in Section |4Jl 
With 
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sig,rms ~ V sou,+ """ ^ sys,+ j^ 
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(A.8) 
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the final result, after shift-and-averaging, is 

It/ >.2 

T'mis.aver =T (T'sig.nnsCv + Av)j + {Tref,rms(y - Av))' 

+ {TZ.ms(y + Av))' + (r-',^,(y - Ay))' 



(A. 10) 



This is roughly a factor of ^/2 lower than for position switching, 
which is simply because for frequency switching the integration 
time on- source is twice as long. 

As in Section lAJl one could again apply a weighting 
scheme, which is more suitable for cases where rsys,off ^ T^^^. 



A.3. Appiication to the example simulations 

To demonstrate that our calibration methods produce the ex- 
pected noise levels, we plot in Fig. lA.ll rnns,aver as calculated 
with Eq. (IA.3I) and dA.lQb using the specific numbers for the 
four temperature phases (which defines the noise in each phase). 
Comparing the theoretical values with the measured RMS levels 
(at two spectral positions) in the figures in Section[3] and IH we 
find that they are consistent. 



Appendix B: Examples of the measurement of Teai 

B.1. Using hot-cold calibration 

Here, we present measurements using the 6. 5 -mm receiver at 
Eff'elsberg. This test was undertaken after a hot-cold measure- 
ment using a continuum backend, such that the feed had already 
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Fig. A.l. Theoretical noise levels for the example simulations. 
The upper panel shows the resulting RMS for position switch- 
ing, and the lower panel contains the frequency switching re- 
sults, which are about a factor of ^f2 lower because of the longer 
eff'ective integration time on- source. 



been exposed to the cold load for some time (about 20 min). The 
receiver was tuned to 43.1 GHz and a bandwidth of 100 MHz 
was chosen. Each load was integrated for 10 min. Furthermore, 
a measurement of the sky was performed afterwards (toward 
zenith). 

Figure IB. II (top panel) shows the integrated raw spectra 
for the three measurements, separated into the cal and non-cal 
phases. The curves reveal the typical footprint of the IF filter 
of the receiving system. Using Eq. (1741) with Thot = 300 K and 
T'coid = 77 K, we computed Giy), which is averaged over the cal 
and non-cal phases (see Fig. IB. II lower panel). 

Dividing the input, P, by the inferred gain curve, G, one ob- 
tains the calibrated data as plotted in Fig. lB.2l (top panel). 

The value of Tcai can be computed via Eq. (l75l) (see Fig. lB.2l 
lower panel). As can be seen, the Tcai spectrum obtained from 
the cold load measurement is much noisier and has a diff'erent 
mean value. This is a very strong indication that the receiver was 
already drifting, probably because the feed was exposed to the 
nitrogen for too long time. Nevertheless, the results illustrate the 
frequency-dependent behaviour of the calibration signal, having 
no 'ripples' or similar defects but a substantial slope considering 
the rather small bandwidth observed. We note that the eff'ective 
bandwidth available is smaller than the nominal 100 MHz owing 
to edge eff'ects. The plot was computed for the inner 50 MHz of 
the spectrum. 



23 



B. Winkel et al.: Unbiased flux calibration methods for spectral-line radio observations 




0.01 0.02 0.03 0.04 



0.05 0.06 
[GHz] 



0.07 



0.08 0.09 
+4.305x10^ 



Fig.B.l. Upper panel: Raw spectra measured during a hot- 
cold experiment with the 6. 5 -mm receiver in spectroscopy mode. 
Each spectrum is the average of a lOmin observation. The sky 
position was observed afterwards by pointing the telescope to 
the zenith. Lower panel: Using the raw spectra measured for 
hot and cold load, one can infer the (calibrated) gain-curve, i.e. 
the transfer function. Residual noise is present but the overall 
shape and even several small features are clearly visible. 

Using the calibrated spectra for the hot-, cold-, and sky- 
measurements, some interesting conclusions can be drawn. As 
for the hot and cold spectra we can assume that 



T-'h,c T-- 

^ sys ~ ^ r? 



+ 7^1 



h,c. 



(B.l) 



one can infer the receiver noise contribution, T^^. From Fig. IB. 21 
(top panel), we have at v = 43.1GHz a Tt. ~ 419 K and 



400 



350 



^300 
250 



200 



150 



n m* ^ M' n w Kw 

M »«iiWHiWi n wii u i|iiiii» m .i 



|> " ** " *|> i w »*iii»* " t'i»>». w »iih*ii<>i M in n iii * Iknt m ^hi It 



Tcold ^hot T'sky 

rpcel rjical T^cal 

^cold ^hot -^sky 



"" *' " *■ '^ " "■'< " * " 'W " ^ »,w..»i M w^N'i. n m ^m m^ n m * ^ * M ¥n < ■■i, a ,., »N *>*fN>,»«,. 



' I '■<*"■* I*" .i...,.ji It ...... , | --*i-f | ff iiiiii -»ini i - i rr uiin.fjvTiiW« 




0.01 



0.02 



0.03 
u [GHz] 



0.04 



0.05 
+4.307x10^ 



Fig. B.2. Upper panel: With the determined transfer function 
(Fig. IB. ll lower panel), one can calibrate the raw spectra (Fig. lB.ll 
upper panel). The individual spectra show a slight slope. Tsky 
is equivalent to the overall system temperature Tsys (which of 
course is a function of observational parameters as telescope 
pointing position, temperature, atmospheric conditions, etc). 
Lower panel: Using the results from the left panel, one can re- 
cover the Teal spectrum. Note that the blue curve T^^}^ is different 
from the two other measurements. (For a discussion of this issue, 
see text.) 



: 0.1 (obtained from an atmospheric 
we calculate r*^^ ^ 263 K . For 



For a z enith opacity of tq = 

model; I Pardoetal.1 l2QQ6l) . 

the sky measurement, the telescope was pointing toward zenith, 

hence AM = 1. However, since weather conditions were very 

good, T might have been somewhat lower, which would lead to 



- sys 



196 K, hence Trx ~ 119K. At 6.5-mm wavelength, the an increase in the inferred Tl^^. We measured an ambient tem- 



system temperature on 'blank' 
approximately by 



sky (in zenith position) is given 



perature of Tamb 
humidity of 57%. 



275 K, air pressure of P = 977 hPa, and 



^ sys 



= T, 



: + T'atm + ^CMB + ^spill,ohrmc • (B.2) 

The antenna gain efficiency can be neglected, as the contribu- 
tions can be treated as constant over the beam pattern. With 
Tsky ~ 148 K and Tcmb ~ 3 K, it follows that Tatm + 7^spiii,ohmic ~ 
26 K. Spill-over can be neglected at this frequency (owing to the 
strong tapering, Tspin < 1 K) and ohmic losses are less than 1 K 
(the reflectivity of the aluminium panels is ~ 0.999). The mea- 
sured value for the Tatm ~ 25 K (at zenith) can be converted into 
the radiation temperature of the atmosphere, r*^^^, using 



^atm = ^ato [1 " exp(-ToAM)] . 



(B.3) 



B.2. Using continuum calibrators 

As for the hot-cold calibration, we present our results of a Tcai 
measurement using continuum calibration sources with accu- 
rately known flux values. Several continuum sources were ob- 
served using the 6-cm receiver in the secondary focus of the 1 00- 
m. The receiving system had an efl'ective bandwidth of 250 MHz. 
The integration time On and Off source was 5min each. For 
the source fluxes of the continuum calibrators, we used new 
values derived by Kraus et al. (in prep.). At 6-cm, one also 
has to correct the spectra for an elevation-dependent antenna 
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Fig. B.3. Measurements of the Tcai spectrum for the 6-cm re- 
ceiver at the 100-m telescope. Two continuum sources (3C 295, 
upper panel; 3C 48, lower panel) were used. Overplotted are val- 
ues obtained from VLB I calibration experiments. Note that the 
right-hand side polarisation signals have been shifted upwards 
by 70 mK for visualisation purposes. 



gairEj. Furthermore, we applied an opacity correction for the 
atmosphere assuming Tzenith = 0.0 16 as a typical value for the 
Effelsberg site at 6-cm. In Fig. IB. 31 we show our results for the 
sources 3C 295 and 3C 48, along with measurements of Tcai from 
VLB I experiments. The two methods lead to consistent results 
with a discrepancy of less than 3%. 



Pohs = P-{ao + aiE\ + a2E\^y 
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